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Abstract 

Purpose The focus of land use modeling in life cycle impact 
assessment has been mainly on taxonomic measures of biodi¬ 
versity, namely species richness (SR). However, increasing 
availability of trait data for species has led to the use of 
functional diversity (FD) as a promising metric to reflect the 
distinctiveness of species; this paper proposes the use of an 
FD index to calculate characterization factors (CFs) for land 
use impacts. Furthermore, we compare the results of the CFs 
to current practice and assess the increase in complexity 
introduced by the use of the new indicator. 

Methods The model proposed is based on data compiled by 
previous regional meta-analysis on SR and FD, in different 
land use types in the Americas. The taxonomic groups included 
were mammals, birds, and plants. Within each study, calculated 
values for FD for different land use types were compared with 
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the natural or close-to-natural state, taken as the reference 
situation. FD values among different land uses were standard¬ 
ized, and CFs were calculated. The final results were then 
analyzed and compared by analysis of variance and post hoc 
tests. A sensitivity analysis was also applied to verify the 
influence on the choice of the reference state. 

Results and discussion The results show that significant dif¬ 
ferences exist between CFs for SR and FD metrics. Across all 
taxa, CFs differ significantly between land use types. The 
results support the use of CF for FD, as a complement to 
current practice. Distinct CFs should be applied for at least six 
groups of land use categories. The choice of reference land use 
type did not significantly alter the results but can be a source 
of variability. A sensitivity analysis evaluating the impact of 
alternate land use types as reference types found only few 
significant changes on the results. 
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Conclusions and recommendations Given the results, we 
believe the use of CFs based on FD can help on the estab¬ 
lishment of possible links between species loss and key 
ecosystem functions, i.e., on the association between the 
midpoint indicator (e.g., biodiversity loss) and the damage 
caused to ecosystem quality, in terms of functions lost. 
Basing CFs on FD is not without challenges. Such indices 
are data hungry (requiring species composition and traits) 
require more complex calculations than current common 
practice, including decisions on the choice of a method to 
calculate FD and the selection of traits. 

Keywords Biodiversity indicator • Functional diversity • 
Global characterization factors • Land use • LCIA • 
Regionalization 

1 Introduction 

The Millennium Ecosystem Assessment (MA 2005) assessed 
the consequences of global ecosystem changes, finding that 
land use change is one of the main direct drivers of biodiversity 
loss, thus impacting the provision of ecosystem services and 
human well-being. According to the Millennium Ecosystem 
Assessment, the conversion rate of natural areas on Earth has 
increased dramatically in the last 50 years especially in tropical 
and subtropical regions, with many biomes have undergone 
up to 50 % change. A direct consequence of such actions 
is the loss of biological diversity, observed by measures 
such as the reduction in population size, increase in homo¬ 
geneity of species composition and increases in species 
extinctions (MA 2005). 

In life cycle impact assessment (LCIA), the focus of 
modeling land use impacts on biodiversity indicators has 
been mainly on taxonomic measures (Curran et al. 2011; 
Vandewalle et al. 2010), such as species richness (SR) 
(Achten et al. 2008; Koellner 2003; Mueller-Wenk 1998; 
Schmidt 2008; Weidema and Lindeijer 2001). However, the 
use of SR alone masks the functional role that species play 
in their habitat (Mouchet et al. 2010). Species are assigned 
an equal weight, regardless of their functional characteristics 
and contributions despite the recognition that species loss 
implies potential changes in the ecosystem related functions. 
For example, carbon storage, nutrient cycling and biotic 
productivity are influenced and affected by activities of 
species. More recently, several researchers have emphasized 
functional diversity (FD) as a more appropriate indicator of 
biodiversity loss in comparison to taxonomic indicators 
because of the association between species traits and eco¬ 
system functions (Diaz and Cabido 2001; Flynn et al. 2009; 
Mouchet et al. 2010; Petchey and Gaston 2006). 

FD is a reflection of the range and value of the quantifi¬ 
able aspects of species—such as feeding behavior, quantity 


of resources consumed, phosphorus uptake, etc.—which are 
measurable at the level of the individual (Petchey et al. 
2009), which if chosen appropriately can act as indicators 
of the roles species play and the rules by which they assem¬ 
ble into communities (Vandewalle et al. 2010) and how they 
influence the way ecosystems operate (Tilman 2001). The 
basis for calculation of FD is a set of functional traits, 
namely the morphological physiological, or behavioral char¬ 
acteristics of organisms (e.g., seed size, leaf area, wood 
density, or plant size, for plants) which can reveal their 
response to environmental pressures (response traits) or the 
effects they have on ecosystem processes (effect traits) 
(Harrington et al. 2010). According to Hooper et al. 
(2002), ecosystem processes are more greatly affected by 
functional differences among species in an ecological com¬ 
munity than their taxonomic richness or composition. 
Theoretically, change in SR does not ensure a change in 
ecosystem processes. The loss of species can be compensat¬ 
ed by other species playing similar functional roles in an 
ecological community. Trait-based measures of FD permit a 
preliminary assessment of redundancy or complementarity 
between species and the functions they play. These relation¬ 
ships cannot be captured by species counting—as currently 
done in LCIA. Any change in species composition has 
consequences on the quantity, range and combination of 
functional traits and therefore has a direct influence on 
the regulation of processes in an ecosystem (Chapin III 
et al. 2000). 

There are two ways to quantify FD: (1) discontinuous 
measures , i.e., group/guild-based, such as the classification 
of species traits according to functional group richness and 
(2) continuous measures , in which species are not divided 
among functional groups (Petchey and Gaston 2006; 
Petchey et al. 2009). Continuous measures all start with 
calculations of the multivariate distance between each pair 
of species in an assemblage. Group-based measures may not 
be appropriate for all ecosystem processes and a large 
amount of decisions and assumptions are required, such as 
where to place the boundaries of each group and the number 
of groups to include. Most importantly, substantial function¬ 
al differences among species are likely to be disregarded in 
some cases when forming these functional groups. The 
second approach consists in calculating the distances between 
species based on their trait values, and then summarizing these 
distances by hierarchical clustering into a dendrogram 
(Petchey and Gaston 2006). Univariate measures based on 
the abundance-weighted range of each trait have also been 
developed (Mason et al. 2003, 2005) after estimating the 
within study and between study variance. 

The calculation of FD index represents a further step on the 
assessment of land use impacts and potentially complement 
current practice in LCIA focusing on SR. In order to increase 
the environmental relevance of the land-use biodiversity impact 
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indicator by better accounting for each species’ role in an 
ecosystem and its stability, the objectives of this study are: 

1. Use a FD index covering several taxonomic levels for 
the calculation of characterization factors for land use 
impacts on biodiversity in LCIA. 

2. Evaluate its influence on biodiversity characterization 
factors compared with current practice. 

3. Assess how far increased complexity and data require¬ 
ments of the new biodiversity indicator are justified by 
improved completeness, reliability, and environmental 
relevance. 


2 Methods 


and land use change (transformation). However, we only 
calculate characterization factors for occupation impacts in 
this model because most of the studies used in the analysis 
did not have information on the previous land use state, which 
is essential for calculating characterization factors for transfor¬ 
mation models. Nor was information on ecosystem recovery 
collected for this study. 

Here, we consider that species not found in a certain land 
use type, in relation to the reference land use, represents 
species loss. The same is true for FD. We emphasize the 
species can be lost without exhibiting changes in FD when 
one or more species playing the same functional role as the 
lost species is retained. 

2.1 Data sources 


The model proposed is based on data compiled by Flynn et 
al. (2009) who analyzed data on SR and FD across land use 
intensification gradients, and studies compiled by Gibson et 
al. (2011). We include three taxonomic groups in this study: 
(1) mammals, (2) birds and (3) plants, mainly due to trait 
data availability. The FD index for each study area was 
completed according to Petchey and Gaston (2002b) using 
presence-absence information for each species. 

Figure 1 shows a schematic view of the data needed for 
the calculation of characterization factors for land use im¬ 
pact assessment based on SR and FD. For methodological 
comparison, we assess differences between the characteri¬ 
zation factors based on FD and SR for different taxonomic 
groups and across land use types. 

Under the framework for LCIA of land use, two types of 
land interventions can be distinguished, land use (occupation) 


The main data sources used in this paper were collected by 
Flynn et al. (2009), in order to carry out an analysis on the 
impact of land use intensity change on SR and FD. In that 
paper, the authors have used three land use categories, ranging 
from natural (e.g., primary forest) to semi-natural (e.g., grazed 
pasture lands) to agricultural (e.g., row crops) systems. 
However, due to the more complex land use classification 
used in LCA, it was essential to investigate the differences 
among more specific land use subcategories. Therefore, addi¬ 
tional studies were added to the original Flynn et al. (2009) 
study to allow the completeness of the analysis carried out in 
this paper. The added studies were cited in Gibson et al. (2011) 
who assessed the impacts on biodiversity in the tropics using a 
meta-analysis of 138 studies. 

Data on functional traits for the additional studies were 
collected from national species databases, such as the Costa 


Fig- 1 Data and methodological 
scheme applied on the calculation 
of characterization factors for 
land use based on SR and FD 
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Rica’s National Biodiversity Institute (INBio), the Portal 
Brasil 500 Passaros , the Animal Diversity Web , and the 
online guide Birds of Surinam. Global databases such as 
BirdLife International , Global Species , NeotropicalBirds , 
and Avian Web have been used when data were not found 
in national databases. 

The mammal meta-analysis included 12 studies conducted 
in Bolivia, Brazil, Pem, Panama, Costa Rica, Mexico, the 
USA, and Canada (Table SI in the Electronic supplementary 
material (ESM). Plant data (Table S2 in the ESM) ranged from 
Costa Rica to eastern Canada. Bird data (Table S3 in the ESM) 
were gathered in Brazil, Costa Rica, Mexico, and the USA. 
Figure 2 presents the spatial location of the data among 
different ecoregions. 

Tables SI, S2, and S3 (ESM) contain more detailed 
information on the data used in the calculations. Data on 
FD, SR, characterization factors for FD and SR, land use 
types and ecoregion are available for each of the studies, for 
all three taxonomic groups. 

2.2 Land use types 

The different land use types present in the studies have been 
aligned with the land use types classification of the 
UNEP/SETAC LULCIA proposal, following Koellner et 
al. (2013). The reference situation chosen to represent the 
potential natural vegetation (PNV) was the most natural or 
close-to-natural state present in each of the studies included 
in the meta-analysis, i.e., primary forest. It is important to 
emphasize that no information on SR of previous natural 
land cover was available in any of the studies and the actual 


land cover was then used. By natural or close-to-natural 
state, we mean land cover types which require a long time 
for restoration (Koellner and Scholz 2007) after disturbance, 
i.e., land use types that represent a stable climax community 
representative of the PNV of the study region given the 
current climatic and edaphic conditions and that show mini¬ 
mal signs of current or recent (i.e., past decades) anthropo¬ 
genic disturbance. For PNV, we use the concept used by Zerbe 
(1998) and Chiarucci et al. (2010) and introduced by Tiixen 
(1956) defined as the hypothetical natural state mature vege¬ 
tation would reach, in a split of second, if all anthropogenic 
interventions would be removed at once. The studies on 
mammals, plants and birds included 13, 13, and 15 different 
land use types, respectively (Table 1). For details, see the 
supplementary Tables SI, S2, and S3 in the ESM. 

2.3 Functional trait data 

The functional traits applied on our analysis were the same 
used in Flynn et al. (2009), because they are the most available 
traits relating species loss with key ecosystem functions 
(Petchey and Gaston 2002a). Accordingly, the main traits 
chosen for this analysis were: (1) resource use and behavioral 
traits for mammals and birds and (2) morphological and 
anatomical traits for plants, related to the species capacity to 
capture resources (Table 2). The same traits information was 
also compiled for the studies cited in Gibson et al. (2011). 

We stress that the objective of this paper is not specifi¬ 
cally to analyze the link between changes in FD and changes 
in ecosystem functions. We recognize that the link between 
species loss and different ecosystem functions needs to be 
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Table 1 Land use and cover classification used in LCA and respective 
codes (Koellner et al. 2013) 


Code 

Land-use-type classification 

[i] 

Forest 

[l.i.i] 

Forest, primary 

[1.1.2] 

Forest, secondary 

[1.2.1] 

Forest, extensive 

[2.2] 

Wetlands, inland 

[3] 

Shrubland 

[4] 

Grassland 

[4.2] 

Pasture/meadow 

[5.1.1] 

Agriculture, arable, fallow 

[5.1.2] 

Agriculture, arable, non-irrigated 

[5.1.2.1] 

Agriculture, arable, non-irrigated, extensive 

[5.1.2.2] 

Agriculture, arable, non-irrigated, intensive 

[5.1.6] 

Field margins/Hedgerow 

[5.2.1] 

Agriculture, permanent crop, non-irrigated 

[6] 

Agriculture, mosaic 

[7.6.3] 

Traffic area, railroad, embankment 

[10.1.2] 

Rivers, artificial 

[10.2.2] 

Lakes, artificial 


further investigated, but this is beyond the scope of the 
current study. 

2.4 FD calculation 

We used Petchey and Gaston’s index of FD in this study 
(Petchey and Gaston 2002b; Petchey et al. 2004, 2009) such 


as in Flynn et al. (2009). Four steps are involved in the 
calculation of Petchey and Gaston’s FD index: (1) construc¬ 
tion of a matrix containing species’ traits values (Fig. 3a); 
(2) calculation of the multivariate distances between species, 
using these trait values (see Fig. 3b); (3) hierarchical clus¬ 
tering of the distance matrix into a dendrogram, which is a 
tree diagram, used to depict the arrangement of clusters (see 
Fig. 3 c); (4) calculation of FD values based on the total 
branch length of the dendrogram, for the species present in a 
particular community. Although abundance-weighted mea¬ 
sures of FD exist, we relied solely on presence-absence data 
in this study due to the lack of species abundance data in 
some of the studies. 

The distance matrix is used to estimate the similarity 
among organisms, in terms of the differences in the trait 
values. We used Gower’s distance metric to calculate the 
distance between species based on their traits, as this dis¬ 
tance metric allows the use of multiple data types, including 
categorical and binary data (Podani and Schmera 2006). 

After obtaining FD values, we normalized values to a 
value of 1 for the reference land use type, in each study 
as in Eqs. (1) and (2) to evaluate FD and SR differences 
within studies. 


II) \ = FD,/FD re f 

(1) 

SR/v = SR,/SR le f 

(2) 


where FD^ and SR N are the normalized values for FD and 
SR, respectively. FD, corresponds to the calculated FD value 
for a certain land use type /, and FD ref is the FD value for the 


Table 2 Functional traits and respective categories used for functional diversity calculation (Flynn et al. 2009) 


Taxonomic group Trait 


Categories 


Birds 

Mass 



Feeding guild 

Carnivore, herbivore, insectivore, and omnivore 


Food type 

Invertebrates, small fruits, seeds, nectar, fish, and generalist 


Foraging location 

Ground, upper canopy, shrub layer, mid-canopy, forage throughout, and aquatic 


Foraging habitat 

Ground, leaves, perch and attack, stems, aerial, water, hover, soar and attack, and other 

Mammals 

Mass 



Feeding guild 

Carnivore, herbivore, insectivore, and omnivore 


Food type 

Invertebrates, fruit, seeds, vertebrates, and vegetation 


Activity 

Diurnal, nocturnal, and either 


Nesting 

Aquatic, arboreal, burrows, multiple, and terrestrial 

Plants 

Leaf area 



Height 



Fruit type 

Meshy and not fleshy 


Fruit length 

— 


Foliage 

Deciduous and evergreen 


Growth form 

Tree, shrub, tall herb, low herb, and grass 


Leguminous 

Legume and not legume 
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Fig. 3 Elements involved on the calculation of Petchey and Gaston’s FD Index: a trait matrix, with trait values v corresponding to species and 
trait t; b distance matrix, with the pair-wise distances between species; and c dendrogram, from which the FD values are calculated 


reference land use type. SR/ and SR ref are the corresponding 
values for SR. 

Contrary to data transformations used by Flynn et al. 
(2009) to compare across studies within land use types, we 
preferred normalization to the computation of standardiza¬ 
tion. In such a way, reference land use types receive a value 
equal to 1, which is more suitable for the calculation of 
occupation impact values. 

We use the natural or most close-to-natural land use type 
as the chosen reference for the normalization procedure 
(Tables S 1 (for mammals), S2 (for plants), and S3 (for birds) 
in the ESM). However, when more than one natural or 
close-to-natural land use types were available within one 
study, we used a sensitivity analysis to verify if a change in 
the results would be observed, by changing the reference 
land use type. 

2.5 Calculation of characterization factors 

We calculated characterization factors based on the natural 
logarithm transformation of normalized FD and SR values to 
adjust the variables to a normal distribution since differences 
within and among studies exist (i.e., area size, sampling 
techniques, etc.). The formulas used for FD and SR are shown 
in Eqs. (3) and (4), respectively. 


CFpp — —Ini FD V ) 

(3) 

CFsr = -M SR n ) 

(4) 


for which FD^ and SR N are the normalized values for FD 
and SR, respectively (Tables SI (for mammals), Table S2 
(for plants), and S3 (for birds) in the ESM). The minus sign 
on Eqs. (3) and (4) was inserted in order to be consistent 
with characterization factors calculated in other models (and 
LCIA impact categories), where impact is expressed in 
positive values. 

An illustration on how the characterization factor 
would be applied in LCIA, to calculate the occupation 
impact (01) is represented in Eqs. (5) and (6), following 
Eqs. (3) and (4), respectively. 


OIfd — CFfd x A occ x t 0 cc (5) 

OIsr = CFsr x A occ x to cc (6) 

A occ represents the area of occupation and t occ the time of 
occupation. 

FD is a complement of SR, i.e., the FD index indicates 
how FD changes (decline or increase) and at which rate, 
when land use change takes place. This can indicate low or 
high redundancy in communities. However, they should not 
be aggregated into one index, once FD index already takes 
into account SR data, but they are rather alternatives to 
represent impacts of land use and land use change on 
biodiversity. 

2.6 Statistical analysis 

Differences between the characterization factors (CF fd x 
CF sr ) were assessed (1) for each taxonomic group and (2) 
among the different land use types existing in the studies. For 
comparison (1), pairwise /Tests and nonparametric Wilcoxon’s 
matched pairs test were applied (Table S4 in the ESM); the 
latter to assess whether deviation from normality affected the 
results. Comparison (2) implied an analysis of variance for all 
taxonomic groups using a significance value of 5 % (0.05). A 
post hoc test was applied to check which groups of data were 
particularly different from each other. We chose Fisher’s least 
significant difference (LSD) test for post hoc comparisons 
among groups and used Mann-Whitney U tests to verify the 
results. We used Statistica 9.0 (StatSoft Inc 2009) to perform 
the statistical analyses. The nonparametric Mann-Whitney U 
test was also applied to verify the results (Table S5 in the 
ESM). Furthermore, we used a weighted analysis, using the 
SR and FD variance as a weight on the previously performed 
analyses (i.e., Mann-Whitney U and LSD tests), to test for the 
influence of differences in sample size and sample variability 
among the land use types. The within-study variance and the 
pooled variance, between studies, was calculated and the 
within and between study variation was compared. The 
weighting procedure was conducted, in order to correct large 
differences between studies variance. For the weighted- 
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variance analyses, the variance of FD and SR values was 
calculated for each of the studies. All the statistical tests were 
then re-run having these values as weight for each data point. 

3 Results 

3.1 Characterization factors for land use impacts 
on biodiversity 

The references used and the calculated values of character¬ 
ization factors for different land use types can be seen in 
Tables SI, S2, and S3 (ESM). These tables present the raw 
values for SR and calculated FD, the respective normalized 
values (SRtv and FD^) and the calculated characterization 
factors based on SR data (CF SR ) and FD (CF fd ). 

3.2 Differences between characterization factors CF fd 
and CF sr for each taxonomic group 

Assessing differences between CF SR and CF fd for all three 
taxonomic groups (without considering different land use 
types) the paired t test revealed a significant difference for 
birds but not for mammals and plants (Table S4 in the 
ESM). 

Similar results were obtained for Wilcoxon’s matched 
pairs test. When this test was applied, taking into account 
different land use types, significant differences were found 
between CF fd and CF SR for “field margins/hedgerows” for 
“birds” and for “agriculture, mosaic” for “mammals.” When 
the data of all groups was used on the analysis, significant 
differences were found for “agriculture, mosaic.” Therefore, 
for these land use types, the results showed that the appli¬ 
cation of different indicators (SR or FD) could lead to 
different conclusions. 

3.3 Difference between characterization factors for different 
land use types 

The analysis of variance showed a degree of significance 
among the CF fd and CF SR values for different land use 
types. A significant difference among the population means 
was found for both CF fd and CF SR , at p levels of 0.0057 
and 0.0003, respectively (Table S6 in the ESM). 

Significant differences were also found among the CFs 
for different land use types, in the results of the post hoc 
Fisher LSD test and the nonparametric Mann-Whitney U 
tests (Table S5 in the ESM). The results show some distinct 
groups, mainly: (1) natural and close-to-natural forest land 
use types, such as “forest, primary,” “forest, secondary,” 
“forest, extensive,” and “forest, intensive” as terrestrial; (2) 
“lakes, artificial” and “rivers, artificial” as water bodies; (3) 
“shrubland” and “field margins/hedgerows”; and (4) 


agriculture land use types, such as, “agriculture, arable, 
nonirrigated, intensive,” “agriculture, mosaic,” “agriculture, 
permanent crop, non-irrigated, intensive,” “pasture/meadow,” 
and “grassland.” Other land use types such as “agriculture, 
arable, non-irrigated,” “agriculture, arable, non-irrigated, 
extensive,” and “agriculture, arable, fallow” did not present 
significant differences towards other land use types. No sig¬ 
nificant differences were observed within land use categories, 
such as forest and agriculture (e.g., between CFs for “forest 
primary” and “forest secondary”). 

The results show that managed landscapes presented a 
detrimental impact in terms of SR and FD loss (Fig. 4). Land 
use types with positive characterization factors (bigger impact, 
in relation to the reference land use types) for both SR and FD 
were: “agriculture, arable, fallow,” “agriculture, arable, non- 
irrigated, intensive,” “agriculture, mosaic,” “agriculture, per¬ 
manent crop, non-irrigated,” “grassland,” “pasture/meadow” 
and “wetlands, inland.” Land use types with negative charac¬ 
terization factors (reduction in SR and FD impact, in relation 
to the reference) were “forest, primary,” “forest, extensive,” 
and “rivers, artificial.” 

The land use type “forest” had unexpected positive CF 
values. “Forest, primary” was taken as reference situation in 
13 studies, “forest, secondary” and “forest” in five, each, and 
other land use types (“shrubland,” “field margins/hedgerows,” 
“grassland,” and “forest, extensive”) in the remaining five 
ones. CFs for “forest, secondary” mainly represents impact, 
when “forest, primary” is taken as reference land use type (see 
Daily et al. (2001) study on Table S3 in the ESM). For “forest, 
extensive”, for which “forest, secondary” was mainly taken as 
reference, both CF SR and CF fd results in values below zero 
(richer in SR and FD). For “forest, intensive,” CF SR values 
were below zero and CF fd mainly above zero. 

For “agriculture, arable, non-irrigated” and “agriculture, 
arable, non-irrigated, extensive,” CF fd values represented 
impact, following the results for other agricultural land use 
types, while CF SR values were mainly below zero. 

CFs for “Traffic area, railroad, embankment” presents 
impacts (positive CFs) for SR and only one negative value 
for CF fd , when “forest, primary” was taken as a reference 
land use type. In that study (Mayfield et al. (2006), the same 
is observed for “pasture/meadow” CF fd and CF SR values: 
they display positive impact, in comparison with “forest, 
primary.” 

“Forest” is taken as reference situation in five studies 
(i.e., CF=0 and therefore, no impact): (1) one bird study 
(Best et al. 1995), (2) two mammal studies (Husband et al. 
2009; Suzan et al. 2008), and (3) two plant studies 
(Sanchez-Merlos et al. 2005a, b). In other four studies, it 
presents positive CFs when compared with the reference 
situations “forest, primary,” “forest, extensive,” and “shrub¬ 
land.” In order to verify the influence of these choices and 
the consequent differences between CF fd and CF SR for each 
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Fig. 4 Box-whisker plot of 
CF sr and CF fd grouped by 
land use types for all taxonomic 
groups (mammals, birds, and 
plants) 
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of the taxonomic groups and for all the data, an analysis of 
sensitivity was carried out. 

3.4 Sensitivity analysis for the choice of reference land 
use types 

Table 3 shows the new reference land use types, in compari¬ 
son to the previous ones, used in the sensitivity analysis. For 
each taxonomic group, few significant changes in the results 
were observed when a different land use type was used as 
reference. The results for the t-tests are displayed in Table S8 
in the ESM. It can be seen that significant differences between 
the means of the characterization factors for FD and SR have 
been found only for “plants.” The results mainly varied due to 
the increase in values of FD and SR for the new reference type 
used for the Mayfield et al. (2006) study, in relation to the 
previous reference used, i.e., most values had a shift from 
negative to positive. For birds, the new reference values 
presented lower SR and FD values for both studies, especially 
in Mas et al. (2004) study, in which most values shifted from 
positive to negative. For different land use types the LSD test 


revealed differences between the CF fd and CF SR before and 
after the sensitivity analysis (Tables S5 and S7 in the ESM). 
The results of the weighted analysis showed quite different 
results, pointing out bigger differences among different land 
use types. For example, differences between crops, pastures, 
and traffic areas could be noted. It seems to be that un¬ 
weighted it is more conservative estimate, because less signif¬ 
icant results appear. 

4 Discussion 

Differences in the magnitude of change in SR and FD between 
broad land use categories, such as forest and agriculture are 
pronounced and therefore characterization factors should be 
differentiated. These results suggest that the integration of the 
ecosystem functions impacted by species loss complements 
existing approaches focusing on SR. 

The resulting characterization factors’ values were as 
expected for both FD and SR indicators: they were higher 
for managed land use types, such as “agriculture, mosaic” 


Table 3 New reference LU types applied in the sensitivity analysis, in comparison to the previous ones, for six studies, including all three 
taxonomic groups 


Taxonomic group 

Study 

Previous reference LU type 

New reference LU type 

Mammals 

Husband (2009) 

Forest ( data 2005) 

Forest (data 2006) 


Sullivan and Sullivan (2006) 

Forest, primary 

Shrubland 

Plants 

Mayfield et al. (2006) 

Forest, primary {forested riverbank) 

Forest, primary (tree-fall gap in forested area) 


Sanchez Merlos et al. (2005b) 

Forest (dry forest ) 

Forest, secondary 

Birds 

Best et al. (1995) 

Floodplain forest 

Grassland ( Upland forest) 


Mas et al. (2004) 

Forest, primary (Forest Beleri) 

Forest, primary (Forest Jrlanda) 
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and lower or negative for close-to-natural or natural 
land use type. However, both indicators did not behave 
equally: this is due to the nonlinear differences on the 
impact caused by biodiversity loss, especially for birds 
and plants. At least for birds, this behavior has also 
been observed in other studies, such as in Flynn et al. 
(2009). Even though, significant differences between 
CF fd and CF SR have been found for only a few land 
use types, FD indicators represent a closer mechanistic 
link between biodiversity loss and resulting impacts in 
ecosystem functioning and stability. 

However, we did not find a significant difference be¬ 
tween sub-land use categories within the same category, 
such as “forest, primary” and “forest, secondary.” The same 
was observed on the analysis done by de Baan et al. (2013). 
The comparison of the values of the FD based CFs to those 
based on SR was a focal point of this study. It should 
however be noted that this only indicates whether or not 
the new approach would lead to different results compared 
with current LCIA practice, while not providing insights 
into which approach is more justifiable and representative. 
It can therefore not be the sole selection criteria between 
both approaches. 

4.1 Classification of land use types 

The first consideration in interpreting these results is the 
classification of land use types to specific categories in each 
study (see Table 1). This involves certain subjectivity on the 
associations made and relies on the description of land use 
types by the authors of the studies. In some of them, only 
few details were given on the specific land covers. The 
category “forest,” for example, is used, when not enough 
information is obtained for more detailed classifications 
such as “forest, primary” and requires only that “an area 
should be covered more than 15 % by trees.” Differences 
between “forest, primary”—“forests minimally disturbed by 
human impact, where flora and fauna species abundance are 
near pristine” (Koellner et al. 2013)— and “forest, 
secondary”—“areas originally covered with forest or wood¬ 
lands, where vegetation has been removed, forest is re¬ 
growing and is no longer in use” (Koellner et al. 2013) — 
have not always been clear in the studies. For other forest 
types “forest” and “forest, extensive,” few data points have 
been found within all data. The differentiation among them 
and other forest types relied on information gathered in the 
studies and was shown to be small. 

4.2 Sample sizes 

The second consideration for these results is the relatively 
small sample sizes for the number of communities in each 
land use type. Ecosystem processes result from the 


interactions among species and the higher the sample size, 
the better it is to define the difference between the simple 
presence or absence of a species versus the functional 
change in ecosystem characteristics (Chapin III et al. 
2000). Small sample sizes might not be adequate to verify 
the variation, especially to what concerns SRx FD measures. 
The analysis of land use types has been conducted across 
different taxa and ecoregions and this could have influenced 
the results. 

4.3 Functional traits 

We used a single set of functional traits for the calculation of 
the FD values. Those traits have been considered by Flynn 
et al. (2009) and suggested to relate to resource capture 
which is the foundation of many ecosystem functions. This 
helps to predict the differences in results between loss of FD 
and SR. 

However, in order to better evaluate the magnitude of 
change due to these functional characteristics, different data 
would be needed. For example, in mammals the trait of 
“offspring size for each species” could drastically change 
as SR declines. However, linking such changes in trait 
diversity to ecosystem functions is challenging, due to rel¬ 
atively few data available within given land use types and 
limitations of such trait data. 

4.4 FD and SR indices 

The results also show the importance of considering com¬ 
plementary information on species presence-absence and 
their functional role in the ecosystem. The quantification 
of aspects of species allows a better detection of specific 
community response to environmental change and impacts 
to ecosystem processes. However, a loss of species may also 
signalize a potential change in the ecosystem functional 
balance, if changes in FD and SR are similar. 

5 Conclusions and recommendations 

Currently, biodiversity indicators used are mainly based on 
taxonomic measures, such as SR. Furthermore, they do not 
take into account important factors such as species ende¬ 
mism, abundance or species vulnerability status, such as 
identified in the IUCN Red List (IUCN 2001) categories 
classification (Souza 2010). These might be crucial for the 
definition of species loss in a local, regional, or global scale 
and also on the capacity of a certain species to recover in a 
certain environment. Although, FD does not necessarily 
include data on these factors, it represents an improvement 
towards a more mechanistic linking between species loss 
and ecosystem processes and stability. 
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The use of FD characterization factors is justified as it 
further allows a more precise detection of community response 
to environmental change and impacts to ecosystem processes, 
although less sensitive to changes and more conservative than 
SR factors. Clearer results would most likely be obtained with 
a larger sample size of communities within each land use type. 
The provision of characterization factors for more detailed 
land use types in LCIA, such as forest types and agriculture 
types remains an issue as the level of detail in the definition of 
different land use categories for LCIA does not correspond to 
biodiversity data availability. Clearly, more data collection 
is needed, also in order to enable the provision of charac¬ 
terization factors for a longer list of different land use 
types. Once the issues with data gaps can be solved, FD 
characterization factors represent a viable and functional 
improvement to current practice. 

Following the recommendations and conclusions of the 
UNEP-SETAC Life Cycle Initiative working group on 
LCIA of land use (sub-group of the working group on 
LCIA of natural resources and land use), the modeling of 
midpoint indicators such as Ecosystem Services Damage 
Potential and Biodiversity Damage Potential was not yet 
made operable. This is mainly due to the need to better 
understand the interactions among species and with their 
environment and to broaden the data compilation on differ¬ 
ent land use types. Therefore, although the current charac¬ 
terization factors are proposed at the midpoint level, it helps 
pointing out possible links between species loss and key 
ecosystem functions. 

The integration of FD in the characterization factors for 
land use impact assessment is not an easy task. Species 
composition needs to be assessed for different land use type 
categories, and data on specific functional traits need to be 
compiled from the literature or directly measured. The final 
characterization factors also involve the calculation of a FD 
index and the existence of different metrics might bring 
differences in the results. The choice of the set of functional 
traits to be used in the calculations may also be a topic of 
concern. However, there is an ongoing effort in order to 
operationalize characterization factors for land use im¬ 
pacts on biodiversity and ecosystem services under the 
UNEP/SETAC Life Cycle Initiative. In line with the re¬ 
sults found, we propose the use of SR indicators, along¬ 
side with FD indices in the assessment of ecosystem 
quality impacts, once they should not be combined in 
one index. This will also allow a better understanding of 
changes in ecosystem functions due to species loss. FD 
based characterization factors allow a clearer link between 
species loss and ecosystem quality indicators, allowing a 
better link between midpoint indicators (biodiversity loss) 
and endpoint indicators (loss of ecosystem functions). 

Beyond biodiversity, land use impact indicators in cur¬ 
rent LCA are covering a growing range of impact pathways, 


such as biotic productivity, carbon sequestration, albedo, 
erosion regulation, or water capturing and filtering/cleaning 
for example. When discussing FD indicators for biodiversi¬ 
ty, it is essential to clearly delimit these service functions 
partially provided by ecosystems from the functional traits 
covered by a biodiversity indicator in order to avoid double 
counting. The functional traits considered in this study are 
resource capture and overall life history strategy, which are 
both related to maintaining the ecosystem (and therefore 
biodiversity). 

Ultimately aiming to increase the reliability (i.e., both 
accuracy and precision) and environmental relevance of 
land-use characterization factors by modeling relevant 
mechanistic relationships between land-use and biodiversity 
impacts this study explored the application of a FD based 
indicator as a complement to SR. Indeed, biodiversity can 
be assessed based on different measures (i.e., taxonomic and 
functional) and a possible recommendation would be the 
possibility to include multi-indicator, capturing many as¬ 
pects of biodiversity and allow the assessment of a range 
of functions. However, this option is not only data demand¬ 
ing, but also requires a better understanding of species in¬ 
teractions in an ecosystem. 

In LCA application, the FD characterization factors can be 
used in the exact same way as the current SR-based factors 
and do not require life cycle inventory data different from the 
latter. The more abundant data availability for plant-SR leads 
to a broad application of this indicator in current LCA prac¬ 
tice, assuming good representativeness for overall biodiversi¬ 
ty. Nevertheless, next to better representing ecosystem 
services and stability, another important advantage of the FD 
indicator is the inclusion of mammals and birds and hence a 
more complete representation of terrestrial ecosystems than 
just plants. However, this strength comes with a price as 
operationalization of this approach for global application in 
LCA is currently limited by the availability of input data 
representing the entire globe and all relevant land use types. 
We thus observe a trade-off between data availability/global 
coverage and representativeness/mechanistic linking to eco¬ 
system functions when concluding on which indicator to 
prefer. Assuming that data availability can be increased 
with further research efforts, FD is a promising biodiver¬ 
sity indicator, which should not be neglected or abandoned 
in future research since this study observed several significant 
differences to SR. 

Acknowledgments We would like to thank the LC Impact (Life Cycle 
Impact Assessment Methods for Improved Sustainability Characterization 
of Technologies, Grant Agreement N.243827, funded by the European 
Commission under the 7th Framework Programme) and SoilTrEC (Soil 
Transformations in European Catchments, Grant Agreement N. 244118, 
funded by the European Commission under the 7th Framework 
Programme) Projects for the financial support provided for the develop¬ 
ment of this research. Further, we are also thankful to M. Gogol-Prokurat, 


4^ Springer 



Int J Life Cycle Assess (2013) 18:1231-1242 


1241 


T. Nogeire, N. Molinari, B.T. Richers, B.B. Lin, N. Simpson, and M.M. 
Mayfield, who through Fabrice DeClerck and Dan Flynn have made part 
of the metadata available for this study 


References 

Achten WMJ, Vandenbempt P, Lemaitre P, Mathijs E, Muys B (2008) 
Proposing a life cycle land use impact calculation methodology, 
6th International Conference on LCA in the Agri-Food Sector, 
Zurich 

Best LB, Freemark KE, Dinsmore JJ, Camp M (1995) A review and 
synthesis of habitat use by breeding birds in agricultural land¬ 
scapes of Iowa. Am Midi Nat 134(1): 1-29 
Chapin FS III, Zavaleta ES, Eviner VT, Naylor RL, Vitousek PM, 
Reynolds HL, Hooper DU, Lavorel S, Sala OE, Hobbie SE, Mack 
MC, Diaz S (2000) Consequences of changing biodiversity. Nature 
405(6783):23A-242 

Chiarucci A, Araujo MB, Decocq G, Beierkuhnlein C, Femandez- 
Palacios JM (2010) The concept of potential natural vegetation: 
an epitaph? J Veg Sci 21:1172-1178 
Curran M, de Baan L, Schryver AMD, Zelm RV, Hellweg S, Koellner 
T et al (2011) Toward meaningful end points of biodiversity in life 
cycle assessment. Environ Sci Technol 45:70-79 
Daily GC, Ehrlich PR, Sanchez-Azofeifa GA (2001) Countryside 
biogeography: use of human-dominated habitats by the Avifauna 
of Southern Costa Rica. Ecol Appl 11(1): 1—13 
De Baan L, Alkemade R, Koellner T (2013) Land use impacts on 
biodiversity in LCA: a global approach. Int J Life Cycle Assess. 
doi:10.1007/sl 1367-012-0412-0 (this issue) 

Diaz S, Cabido M (2001) Vive la difference: plant functional diversity 
matters to ecosystem processes. Trends Ecol Evol 16(ll):646-655 
Flynn DFB, Gogol-Prokurat M, Nogeire T, Molinari N, Richers BT, 
Lin BB, Simpson N, Mayfield MM, DeClerk F (2009) Loss of 
functional diversity under land use intensification across multiple 
taxa. Ecol Lett 12(1 ):22—33 

Gibson L, Lee TM, Koh LP, Brook BW, Gardner TA, Barlow J, Peres 
CA, Bradshaw CJA, Laurance WF, Lovejoy TE, Sodhi NS (2011) 
Primary forests are irreplaceable for sustaining tropical biodiver¬ 
sity. Nature 478:378-383 

Harrington R, Anton C, Dawson TP, Fd B, Feld CK, Haslett JR, 
Kluvankova-Oravska T, Kontogianni A, Lavorel S, Luck GW, 
Rounsevell MDA, Samways MJ, Settele J, Skourtos M, 
Spangenberg JH, Vandewalle M, Zobel M, Harrison PA (2010) 
Ecosystem services and biodiversity conservation: concepts and 
glossary. Biodivers Conserv 19(10):2773-2790 
Hooper DU, Solan M, Symstad A, Diaz S, Gessner MO, Buchmann N, 
Degrange V, Grime P, Hulot F, Mermillod-Blondin F, Roy J, 
Spehn E, van Peer L (2002) Species diversity, functional diversity, 
and ecosystem functioning. In: Loreau M, Naeem S, Inchausti P 
(eds) Biodiversity and ecosystem functioning: synthesis and per¬ 
spectives. Oxford University Press, London, pp 195-208 
Husband TP, Abedon DH, Donelan E, Paton P (2009) Do coffee 
dominated landscapes support mammal biodiversity? In: World 
Agroforestry Centre (Hrsg.), 2nd World Congress of Agroforest¬ 
ry, Agroforestry - The Future of Global Land Use. Book of 
Abstracts, 2nd World Congress of Agroforestry. World Agrofor¬ 
estry Centre, Nairobi, Kenya 

IUCN (2001) IUCN Red List Categories and Criteria. Version 3.1. 

IUCN Species Survival Commission, Gland 
Koellner T (2003) Land Use in Product Life Cycles and Ecosystem 
Quality. Peter Lang, Frankfurt am Main, 271 pp 
Koellner T, de Baan L, Brandao M, Mila i Canals L, Civit B, Goedkoop 
M, Margni M, Weidema BP, Wittstock B, Mueller-Wenk R (2013) 


Principles for life cycle inventories of land use on a global scale. Int J 
Life Cycle Assess. doi:l 0.1007/s 11367-012-0392-0 
Koellner T, Scholz RW (2007) Assessment of land use impacts on the 
natural environment. Part 1: An analytical framework for pure 
land occupation and land use change. Int J Life Cycle Assess 
12:16-23 

MA (2005) Ecosystems and Human Well-being: Biodiversity Synthe¬ 
sis. World Resources Institute, Washington, DC 
Mas AH, Dietsch TV (2004) Linking shade coffee certification to 
biodiversity conservation: butterflies and birds in Chiapas, Mex¬ 
ico. Ecol Appl 14(3):642-654 

Mason NWH, MacGillivray K, Steel JB, Wilson JB (2003) An index of 
functional diversity. J Veg Sci 14(4):571-578 
Mason NWH, Mouillot D, Lee WG, Wilson JB (2005) Functional 
richness, functional evenness and functional divergence: the pri¬ 
mary components of functional diversity. Oikos 111(1): 112—118 
Mayfield MM, Ackerly D, Daily GC (2006) The diversity and conser¬ 
vation of plant reproductive and dispersal functional traits in 
human-dominated tropical landscapes. J Ecol 94(3):522-536 
Mouchet MA, Villeger S, Mason NWH, Mouillot D (2010) Functional 
diversity measures: an overview of their redundancy and their 
ability to discriminate community assembly rules. Funct Ecol 
24(4): 867-876 

Mueller-Wenk R (1998) Land use the main threat to species: How to 
include land use in LCA. 64. Universitat St. Gallen, St. Gallen 
Petchey OL, Gaston KJ (2002a) Extinction and the loss of functional 
diversity. Proc R Soc Lond B 269:1721-1727 
Petchey OL, Gaston KJ (2002b) Functional diversity (FD), species 
richness and community composition. Ecol Lett 5(3):402^411 
Petchey OL, Hector A, Gaston KJ (2004) How do different measures 
of functional diversity perform? Ecology 85(3):847-857 
Petchey OL, Gaston KJ (2006) Functional diversity: back to basics and 
looking forward. Ecol Lett 9(6):741-758 
Petchey OL, O’Gorman EJ, Flynn DFB (2009) A functional guide to 
functional diversity measures. In: Naeem S, Bunker DE, Hector 
A, Loreau M, Perrings C (eds) Biodiversity, ecosystem function¬ 
ing, and human wellbeing: an ecological and economic perspec¬ 
tive. Oxford University Press, Oxford, p 384 
Podani J, Schmera D (2006) On dendrogram based measures of func¬ 
tional diversity. Oikos 115(1):179—185 
Sanchez-Merlos D, Harvey CA, Grijalva A, Medina A, Vilchez S, 
Hernandez B (2005 a) Diversidad, composicion y estructura de 
la vegetacion en un agropaisaje ganadero en Matiguas, Nicaragua. 
Rev Biol Trop 53(3-4):387-414 

Sanchez-Merlos D, Harvey CA, Grijalva A, Medina A, Vilchez S, 
Hernandez B (2005b) Caracterizacion de la diversidad, densidad 
y estructura de la vegetacion en un paisaje fragmentado de bosque 
seco en Rivas, Nicaragua. Rev Rec Nat Ambiente 45:91-104 
Schmidt JH (2008) Development of LCIA characterization factors for 
land use impacts on biodiversity. J Cleaner Prod 16(18): 1929^42 
Souza DM (2010) Proposta de um modelo de caracterizagao de 
impactos do uso da terra, segundo indicadores de biodiversidade, 
em AICV: calculo de fatores de caracterizagao para ecorregides 
brasileiras. Ph.D. thesis, Universidade Federal de Santa Catarina, 
Florianopolis, 309 pp 

StatSoft Inc (2009) Statistica: data analysis software system. StatSoft. 
Inc, Tulsa 

Sullivan TP, Sullivan DS (2006) Plant and small mammal diversity in 
orchard versus non-crop habitats. Agr Ecosyst Environ 116(3^4): 
235-243 

Suzan G, Armien A, Mills JN, Marce E, Ceballos G, Avila M, Salazar- 
Bravo J, Ruedas L, Armien B, Yates TL (2008) Epidemiological 
considerations of rodent community composition in fragmented 
landscapes in panama. J Mammalogy 89(3):684-690 
Tilman D (2001) Functional diversity. In: Levin SA (ed) Encyclopedia 
of biodiversity. Academic, New Jersey 


Springer 



1242 


Int J Life Cycle Assess (2013) 18:1231-1242 


Tiixen R (1956) Die heutige potentielle naturliche Vegetation als Gegenstand 
der Vegetationskartierung. Angew. Pflanzensoziol (13):5-42 
Vandewalle M, de Bello F, Berg M, Bolger T, Doledec S, Dubs F, Feld 
C, Harrington R, Harrison P, Lavorel S, da Silva P, Moretti M, 
Niemela J, Santos P, Sattler T, Sousa J, Sykes M, Vanbergen A, 
Woodcock B (2010) Functional traits as indicators of biodiversity 
response to land use changes across ecosystems and organisms. 
Biodivers Conserv 19(10):2921—2947 


Weidema B, Lindeijer E (2001) Physical impacts of land use in product 
life cycle assessment. Final report of the EURENVIRON- 
LCAGAPS sub-project on land use, Department of Manufactur¬ 
ing Engineering and Management, Technical University of Den¬ 
mark, Lyngby 

Zerbe S (1998) Potential natural vegetation: validity and applicability 
in landscape planning and nature conservation. Appl Veg Sci 
1:165-172 


4^ Springer 



